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The  Effect  of  Moisture  on  Carbon  Fiber  Reinforced  Epoxy  Composites. 
I Diffusion. 


The  observation  by  earlier  investigators  that  the  flexural  and 
shear  strengths  carbon  fiber  reinforced  composites  may  deteriorate 
at  eunbient  storage  condition  made  it  necessary  to  investigate  these 
phenomena  more  closely  in  order  to  get  a better  understanding,  and 
ultimately,  to  be  able  to  predict  how  these  materials  would  behave 
in  longtime  service  and  storage  conditions. 

The  work  described  in  this  report  is  a part  of  this  investi- 
gation. It  describes  the  diffusion  of  moisture  in  carbon  fiber 
composites  which  may  serve  as  a basis  for  prediction  of  longtime 
composite  changes  in  various  environments. 

This  program  was  funded  by  the  Naval  Air  Systems  Command 
(Task  No.  A3200000010123)  during  the  periods  of  1 July  1974  to 
1 January  1976. 
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INTRODUCTION 


The  use  of  carbon  fiber  reinforced  composites  for  primary 
and  secondary  structural  components  in  aircraft  construction 
promises  not  only  benefits  in  weight  savings,  and  with  it,  maneu- 
verability and  range  of  the  aircraft,  but  also  a considerable 
savings  in  labor  cost.  This  savings  in  labor  cost  more  than 
compensates  for  the  higher  cost  of  the  raw  materials. 

Serious  consideration  is  given  today  to  the  replacement  of 
various  components  in  aircraft  with  carbon  fiber  composites.  For 
instance,  the  Lockheed  Aircraft  Corporation  has  decided  to  use 
carbon  fiber-epoxy  composites  in  the  future  for  the  vertical 
stabilizer  of  their  1011  commercial  airliner.  It  is  true  that 
while  the  demands  of  commercial  aircraft  are  less  severe,  than 
those  of  military  aircraft  the  above  use  also  indicates  that 
confidence  in  these  advanced  composites  has  increased  over  the 
past  few  years.  Therefore,  it  becomes  even  more  necessary  to 
learn  as  much  as  possible  about  these  materials,  such  as  their 
behavior  in  various  environments  and  the  change  of  properties 
with  time. 

This  report  describes  part  of  an  investigation  to  highlight 
certain  aspects  of  environmental  property  changes  in  carbon  fiber 
reinforced  composites.  There  are  a number  of  mechanisms  that  can 
lead  to  changes  of  composite  properties.  In  this  investigation 
emphasis  was  given  to  the  effect  of  moisture  on  these  composites, 
since  the  rate  of  moisture  absorption  at  ambient  conditions  is 
considerably  higher  than  the  rate  of  any  other  mechanism  that  may 
lead  to  property  changes  in  benign  environments.  Specifically, 
the  subject  of  this  report  is  the  diffusion  process  of  moisture 
in  high  performance  carbon  fiber  reinforced  composites.  Diffusion 
coefficients  were  determined  on  several  composite  systems  as  a 
function  of  the  relative  humidity  and  temperature.  These  data 
serve  as  basis  for  estimating  the  moisture  distribution  inside 
the  composites  as  functions  of  composite  thickness,  time  and 
temperature.  The  concomitant  composite  property  changes  will  be 
described  in  a subsequent  technical  report. 

BACKGROUND  INFORMATION 

The  effect  of  moisture  on  high  temperature  properties  of 
carbon  fiber  epoxy  composites  was  discovered  jointly  by  investi- 
gators of  the  Fiber ite  and  Aerospace  Corporation  and  is  well 
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documented  in  a summary  report  by  J.  Hertz  [l] . C.  C.  Browning  [2] 
at  the  Air  Force  Materials  Laboratory  (AFML)  compared  composites 
with  different  fiber  orientations  and  showed  that  the  properties 
which  are  most  strongly  affected  by  moisture  are  the  matrix 
controlled  properties. 

In  earlier  repcil-ts,  R.  Simon  et  al.  [3]  at  the  Naval  Surface 
Weapons  Center  had  investigated  the  effect  of  water  on  various  ^ 
carbon  fiber  composites  and  concluded  that  mechanical  property 
changes  are  affected  by  both  the  nature  of  the  carbon  fiber  and 
the  type  of  resin  used. 

The  longterm  degradation  effects  on  fiber  reinforced  composites 
can  be  caused  by  various  mechanisms.  Which  mechanisms  predominate 
depends,  of  course,  on  the  type  of  materials  that  constitute  the 
composite.  In  fiber  reinforced  composites  a degradation  of  the 
resin,  the  fiber  or  the  interfacial  bonding  will  have  a pronounced 
effect  on  various  composite  properties. 

There  are  a number  of  mechanisms  that  can  lead  to  composite 
property  changes  which  may  be  differentiated  as  reversible  and 
irreversible  changes. 

Chemical  and  mechanical  matrix  degradation  such  as  photo- 
oxidation, thermal  degradation  (in  the  presence  or  absence  of 
air),  hydrolytic  polymer  chain  cleavage,  high  energy  radiaiton, 
stress  cracking,  microbial  degradation,  and  thermal  [4]  and 
mechanical  cycling  may  lead  to  irreversible  property  changes  of 
fiber  reinforced  composites. 

We  define  reversible  property  changes  as  those  changes  that 
can  be  reversed  to  the  original  properties  (at  least  within  the 
limitations  of  the  sensitivity  of  the  measurements) . 

A trivial  case  of  a reversible  property  change  is  a change  due 
to  temperature,  as  long  as  the  temperature  change  is  not  too  high 
which  could  lead  to  an  irreversible  degradation. 

Another  reversible  change  in  composite  properties  is  observed 
when  small  molecules  such  as  moisture  or  organic  molecules  penetrate 
the  matrix  and  plasticize  it.  The  effect  of  such  a plasticization 


Ul  J.  Hertz,  Final  Report  NASA8-27435,  June  1973  (NASA  Contract). 

[2J  C.  E.  Browning,  28th  Annual  Conference  of  the  Society  of  the  • 

Plastic  Industry,  Feb  1973  (Washington,  D.C.);  Proceedings  15A. 

[3]  R.  A.  Simon  and  S.  P.  Prosen,  NOLTR  68-132,  Oct  1968; 

M.  L.  Santelli  and  R.  A.  Simon,  NOLTR  70-258,  Feb  1971.  » 

[4]  E.  L.  McKagne,  Jr.,  J.  E.  Halkias  and  J.  D.  Reynolds,  "Moisture 

in  Composites:  "Effects  of  Supersonic  Service  on  Diffusion," 

J.  Comp.  Materials,  9,  2,  (1975). 
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is  a reduction  in  the  glass • transition  temperature  (T^)  of  the 
matrix  and  also  a reduction  in  the  matrix  modulus  over  a wide 
temperature  range.  By  heating  and/or  evacuating  of  such  materials 
one  can  remove  the  moisture  again  and  the  original  strength 
properties  may  be  regained. 

The  fact  that  this  reaction  is  reversible  is  of  little  practical, 
value,  since  the  strength  properties  of  a large  structure  cannot  be 
easily  restored.  On  the  other  hand,  the  composite  strength  is 
expected  to  deteriorate  only  to  a certain  level  depending  on  the 
equilibrium  conditions,  which  are  governed  by  the  environment, 
especially  by  the  relative  hximidity. 

The  Experimental  section  will  briefly  describe  the  materials 
we  have  investigated  and  the  techniques  for  measuring  the  diffusion 
coefficients.  The  following  Discussion  section  will  focus  on  the 
moisture  diffusion  in  composite  materials,  while  the  resulting 
changes  in  mechanical  properties  will  be  described  in  a subsequent 
report. 


EXPERIMENTAL 

Materials 


A.  Resins 


1.  Narmco  5208.  This  resin  was  obtained  from  the  Whittaker 
Corporation,  Narmco  Division.  It  is  a one-component  system  'resin 
plus  curing  agent) . 

2.  DER  332/DADS.  The  resin  system  consists  of  100  parts.of 
DER  332  (a  diglycidyl  ether  of  bisphenol  A)  and  36  parts  of  DADS 
(4,4'  diaminodiphenyl  sulfone) . 

3.  Epon  1031/NMA.  This  resin  consists  of  100  parts  of  Epon 
1031  [l, 1 , 2 , 2-tetra (p-glycidyloxyphenyl)  ethane  ],  77  parts  of 
Nadic  methyl  anhydride,  and  1 part  of  BDMA  (benzyldimethyl  amine) . 

B.  Fibers 


1.  Thornel  300.  This  material  was  obtained  from  Union 
Carbide  Corporation  and  consists  of  a continuous  strand  of  3,000 
filaments  which  have  an  epoxy  sizing  for  better  handleability . 

(For  properties  reported  by  the  manufacturer  see  Table  1.) 

2.  HMS  fibers.  This  material  is  a tow  with  10,000  continuous 
filaments  manufactured  by  the  Hercules  Corporation.  The  fibers 
had  no  sizing.  (For  properties  reported  by  the  manufacturer  see 
Table  1.) 
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The  composites  were  made  from  8 plies  of  9-inch  wide  prepregs, 
which  were  layed  up  unidirectionally  and  molded  (in  a vacuum  bag) 
under  100  psi  pressure.  The  details  of  the  prepregging  operation 
and  the  composite  curing  conditions  will  be  described  in  a forth- 
coming technical  report  that  will  deal  with  the  moisture  and 
temperature  effects  on  the  mechanical  properties  of  these  composites. 


Table  1:  Carbon  Fiber  Properties 


» 


Properties 

Hercules  HMS 

Thornel  800 

Strand  Modulus,  psi  x 10^ 
(G  Pasc. ) 

50  - 55 
(345  - 379) 

33  - 34.5 
(227  - 238) 

Strand  Break  Strength,  psi  x 10^ 
(G  Pasc. ) 

340  (2.34) 

361  (2.49) 

Density,  g/cc  (lot  average) 

1.85  - 1.90 

1.74  - 1.78 

Filament  Diameter,  microns 

7.2  - 7.5 

8 

Number  of  Filaments 

10,000 

3,000 
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Measurements 

A.  Equilibrium  Concentrations.  To  determined  the  moisture 
equilibrium  concentrations  of  the  resins  and  composites  fine 
machine  shavings  were  obtained  and  dried  for  48  hours  at  120“C  in 
vacuum.  The  dried  shaving  samples  were  stored  over  magnesium 
perchlorate  desiccant  before  use.  One  gram  samples  were  prepared 
in  individual  crucibles.  The  sample  exposure  during  the  weighing 
operation  was  kept  below  3 minutes.  The  samples  were  then  exposed 
to  30%  and  80%  RH,  at  30®,  45®,  60®,  and  75°C. 

The  equilibrium  concentration  was  defined  as  the  concentration 
when  the  weight  increase  vs.  time  curve  leveled  out,  which 
generally  occurred  within  a few  days.  The  percent  weight  gain 
after  two  weeks  was  designated  the  equilibrium  concentration. 

B.  Fractional  Weight  Increase.  5 x 5 cm  square  samples  were 
machined  from  the  composite  panels  and  sealed  around  the  edges  with 
a strip  of  aluminum  foil  (adhesively  bonded) . These  specimens  were 
dried  for  72  hours  at  125®C  before  exposure  in  the  humidity 
chambers  (30%  and  80%  RH , at  30®,  45®,  60®  and  75®C) . The  plate 
samples  were  removed  after  various  time  intervals  and  the  fractional 
weight  increase  (%) (%) ] was  plotted  vs.  square  root  of 

time,  where  is  the  percent  weight  increase  after  the  time(t) 
and  is  the  equilibrium  value. 

Humidity  Chambers 

The  samples  were  put  into  a sample  holder  and  placed  in  a 
wide-mouth  screw  cup  container  (23.5  cm  high,  11.5  cm  diameter). 
Various  salt  solutions  (2  cm  high,  with  undissolved  solute)  served 
to  keep  the  relative  humidity  constant.  The  exposure  chambers 
were  closed  and  submerged  into  4 water  baths  held  at  30®,  45®, 

60®,  and  75®C  (+.2®C) . 


DISCUSSION  AND  RESULTS 


A.  Materials 


Considerations  of  aerodynamic  heating  which  predicted  that  parts 
of  the  aircraft  may  reach  350®F  (177®C)  during  supersonic  flight 
made  it  mandatory  that  the  composite  materials  to  be  used  for  these 
aircraft  be  stable  at  this  temperature.  Composites  with  177®C 
performance  capability  require  resin  matrices  with  a glass  transition 
temperature,  Tg,  > I77<>c,  i.e.,  resins  that  have  no  considerable 
loss  of  modulus  below  this  temperature.  Because  absorption  of 
moisture  leads  to  a plasticization  of  the  matrix  it  was  expected 
that  resins  with  low  moisture  affinity  would  show  a better  strength 
retention  than  those  with  high  moisture  saturation  values. 

This  work  was  restricted  to  the  investigation  of  six  composite 
systems  made  from  two  different  carbon  fibers  and  three  resins. 
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The  resins  were:  Narmco  5208,  DER  332/DADS  and  Epon  1031/NMA. 
The  fibers  were  Thornel  300  and  HMS  fiber. 

During  the  original  screening  of  resins  by  torsional  braid 
analysis  for  their  thermomechanical  behavior,  we  found  that  several 
other  resins  also  fulfilled  the  temperature  requirements.  Some  of 
the  properties  of  these  resins  are  listed  in  Table  2. 

Two  of  these  resins,  Fiberite  904  and  X-801,  are  no  longer 
commercially  available.  The  choice  of  resins  for  this  investigation 
was  somewhat  arbitrary.  For  purposes  of  scientific  understanding 
those  resins  whose  chemical  composition  was  known  were  preferred. 
Nevertheless,  Narmco  5208,  a proprietary  commercial  material  was 
also  included  because  of  its  attractive  properties.  (Hercules  3501 
appeared  similarly  attractive;  we  intend  to  investigate  this 
material  at  a later  time.) 

The  choice  of  carbon  fibers  was  based  mainly  on  the  similarity 
in  tensile  properties  and  on  their  dissimilarity  in  modulus  (the 
manufactures'  data  are  given  in  Table  1). 

B.  Effect  of  Moisture  on  Composite  Strength 

The  change  of  composite  strength  as  a function  of  temperature 
and  moisture  will  be  the  subjects  of  two  forthcoming  technical 
reports.  It  will  suffice  here  to  summarize  what  might  be  expected 
if  a composite  is  exposed  to  humidity. 

At  ambient  temperature  the  rate  of  degradation  of  composite 
strength  due  to  oxidation  or  hydrolysis  is  very  small  compared  to 
that  of  moisture  absorption.  Thus  it  should  be  possible  to  predict 
the  ultimate  strength  and  modulus  changes  which  are  reached  when 
the  moisture  content  in  the  composite  is  in  equilibrium  with  its 
surrounding.  This  equilibrium,  of  course,  is  a function  of  the 
temperature  and  of  the  relative  humidity.  The  degradation  of 
composite  properties  is  expected  to  be  a function  of  the  absorbed 
moisture : 


10 


11  Itr^~  ' 


r 


M 


rn 


H-l 

o 

19 

•H 

« 

a 

o 

(U 


CM 

(V 

iH 

Xi 

la 


0) 

» 

A! 

U 

0 

01 


(0 

a 


M 


0) 

p 


NSWC/WOL/TR  76-7 


Eh 

O 

CN 

00 

ro 

ro 

fO 

00 

00 

m 

1 

1 

00 

CN 

CN 

rH 

CN 

VO 

CN 

+ 

rH 

+ 

+ 

+ 

n. 

+ 

+ 

+ 

<N 

ro 

in 

ro 

1 

CN 

in 

o 

CN 

m 

tH 

iH 

00 

1 

m 

rH 

CM 

(N 

CN 

CN 

CN 

CN 

CN 

<N 

CN 

<N 

CN 

rH 

r- 

O 

in 

O 

O 

O 

in 

in 

CN 

rH 

r- 

r- 

VO 

O 

O 

rH 

o 

o 

00 

rH 

rH 

rH 

CN 

CN 

CN 

n. 

CN 

CN 

CN 

rH 

in 

in 

m 

in 

rH 

rH 

rH 

rH 

m 

o 

o 

r- 

o 

m 

rH 

in 

00 

in 

rH 

O 

<0 

rH 

m 

m 

CN 

o 

VO 

rn 

in 

rH 

rH 

cr> 

CN 

CN 

CO 

m 

CN 

CN 

CN 

cn 

CO 

CM 

<0 

tH 

•o 

•o 

at 

c 

■P  M iH 

Id 

C P 

a 

a\ 

00 

rH 

rH 

in 

rH 

CO 

u 

at  4J 

o 

rH 

VO 

CO 

rH 

u 10 

tX 

• 

• 

• 

• 

• 

• 

• 

• 

• 

c 

P -H 

u 

CO 

in 

VO 

r*- 

rH 

CM 

•tH 

0)  0 

•tH 

09 

IX  £ 

a 

0) 

K 

>■ 

r- 

rH 

o 

r-* 

in 

CO 

OV 

>1 

■P 

m 

in 

CO 

in 

o 

CO 

X 

•H 

u 

CN 

CN 

CN 

CN 

CN 

CN 

CO 

CN 

CN 

Q 

10 

u 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o< 

c 

\ 

rH 

rH 

rH 

iH 

rH 

rH 

r'K 

*■  rH 

U 

0) 

o< 

Q 

c 

o 

CO 

in 

CO 

o 

(N 

in 

00 

•o 

a> 

o 

o 

o 

o 

O 

CO 

CN 

iH 

iH 

rH 

tH 

rH 

rH 

rH 

Id 

£ 

+ 

+ 

+ 

+ 

+ 

U 

+ 

u 

u 

u 

u 

o 

u 

u 

u 

u 

o 

o 

o 

CJ  u 

u 

u 

e 

u 

U 

o 

U 

U 

U U 

in 

o 

u 

0 

o 

o 

O 

O 

o 

o o 

0 

o 

o 

o 

o 

o 

0 

o 

o o 

CN 

O 

o 

o 

O 

O 

o 

o 

O 

O O 

CO 

O 

o 

O 

O 

O 

in 

O 

o 

rH 

o 

in 

in 

O 

VO 

rH 

CN 

tH 

m o 

CN 

o 

rH 

O 

in 

O 

CN 

in 

CN  VO 

CM 

o^ 

tH 

tH 

IN 

iH  IN 

IN 

IN 

rH 

CN 

rH 

CN 

rH  CN 

at 

•P 

P 

P 

P 

P 

P 

tH 

•P 

P 

•P 

■P 

P 

Id 

Id 

Id 

p p 

(d 

P 

Id 

P 

P 

P 

P 

P 

P P 

u 

Id 

Id 

Id 

Id 

Id 

Id  Id 

Id 

Id 

Id 

Id 

<0 

Cd 

(d  cd 

>1 

u 

hr 

hr 

hr 

hr 

£ 

hr 

hr 

hr 

hr 

hr 

hr 

hr 

hr 

hr 

hr 

hr 

VO 

VO 

VO 

CN 

m 

VO 

rH 

CN 

rH 

CO  ^ 

rH 

CN 

CO 

CO 

in 

^ CM 

CN 

at 

in 

at 

w 

CN 

p 

CN 

P 

u 

CO 

•H 

in 

■»r 

o 

m 

00 

in 

•H 

l-t 

c 

CO  CO 

p 

rH 

m r*> 

U 

tH 

X 

tH  m 

P 

ON 

D 

•H 

(0 

S 

O VO 

r « 

S 00 

O 

o 

o w 

rH  VO 

U-H 

0) 

PC  Q 

.Q 

00  < 

' Q 

» o 

00 

z 

00  Q 

ia 

m 

«o 

at 

u 3 

•H 

m 

1 sc 

9 < 

CN 

1 

O 

1 < 

•H 

o 

O £ 

am 

K 

o o 

a 

as 

X Z 

S Q 

z in 

X 

ti 

X D 

a 

Ov 

rH  Z 

sm 

11 


NSWC/WOL/TR  76-7 


CO 

Ip 

0 

Tl 

c 

0) 

3 

p 

0 

3 

P 

0) 

P 

O' 

c 

4 

•H 

P 

c 

E 

0) 

0 

4 

a 

•P 

E 

>1 

•a 

0) 

CO 

IQ 

p 

0) 

p 

0) 

0) 

•p4 

c 

p 

P *0 

0) 

0) 

4 

a 

0 £ 

*H 

fH 

£ 

•P 

0) 

0)  4 

p 

>1 

P 

> 

0) 

a 

5 

c 

0) 

p 

X -C 

0) 

CO 

'0 

CO 

0)  84 

£ 

4 

<u 

a 

CO 

03 

0)  — £ 

> 

4 

T3 

3 

c 

£ U 0 

•H 

p 

0) 

fp 

4 

0 -P 

4J 

0) 

C 

3 

0 P 

E 

*P 

»o 

P 

p P 

rH 

ip 

0 

4 

«p  0 

0) 

•O 

0) 

E 

fH 

0) a 

u 

• £ 

•a 

p 

> 4 

C 

3 a p 

dP 

•H 

0) 

-^•p 

<V 

P 0 £ 

in 

Q 0) 

4 0 O' 

C 

3 

(0  CO 

p in  -P 

0 -H 

P 84  CO 

0)  ro  4 

0 *H 

to 

0 

a M 

c 

E -H 

>1^ 

0 

g .p 

<T) 

C C 

0)  £ 

0)  • CO  O' 

< 

4 4 

a 

0) 

C 

P 0)  £ 

0) 

Z -P 

E 

£ 

0 

0)  • 0 -P 

iH 

rH 

T3 

0) 

O'  P 

•r4 

> 0)  -P  £ £ 

C 9>  • C 

4 0)  • £ a V 

U Si  Ouw  ^ 

0)  O4  M O >1 

a (0  0)  o E u x: 

6 -H  C U O 0)  -u 

0)  ja  o Tt  ^4  V 

4J  «M  0)  •'■I  0)  E 

IW  i-t  4J 

e o 3-H  c'04'»-i 

O CO  >>l  H >1  ^ O 
O M 0)  10  j::  -P 

Mfl).H^0)C<H0) 


0) 

p 

c 

p 

p 

0 

4 

0) 

• 

£ 

•H 

0) 

0 

3 

E 

in 

CO 

p 

a 0) 

£ 

a p 

0) 

0 

4 

£ 

4 

E 

p 

•H 

0 

s 

0) 

CO 

E 

0 £ 

0) 

0 

>1 

•P 

4 £ 

p 

CO 

P 

•P 

•0 

£ 

u 

£ 

P 

•o 

> 

0) 

0) 

3 

0 

P 

0) 

0) 

P 

P 

•H 

£ P 

1-4  ip 

•p 

> 

£ 

p 

91 

Vi 

<P 

0) 

•H 

c: 

0 

P 

T5 

P 

•P 

0) 

0) 

Sfl 

E 

C 

4 

p 

iH 

•iH 

0) 

£ 

0 

to 

0 

3 

Q 

a 

0 

0 

0) 

£ 

4 

Ip 

u 

•P 

P 

p <p 

G 

0 

-a 

•P 

•H  -M  C 4 
iH  4 V4 
0)  CO  £ 0) 

£ CO  u a 


s:  >1-H  V44£M£COOQ4 
^■PCfo  S3cn44g 
40)0)  >-11-1  4-1  i-ieo) 

£ >i4  >i>iX 
•O  ^ a£  4J  £ £ -H  £ E 

0)  >|-H  C -P  E 0)  3 0) 
C'O'OTJWa'O  '0£E£ 

•H-r4  00)EE«)O0)-P-HP 

go  c -H  -iH  ^ -iH  X 

r-^^lOr^\r^C0  4C0 


go  c 

>i-H 


0)  iH  E a a)  -H  ao 
4J  tp  4 a a-o  avo 
0)  -H  -H  3 X 4 3 


•H  X 

CO  4 

a**^  E 
a 

3 O' 


£ O E «0 

O'  0)  3 0)  3 o 

•H  P £ E O --C  C 

4 c P -H  -r4  3 

P r4  X > T}  0) 

P Q 'p  4 P O > 

CO  a O E 0)  E 4 

m £ 

4 0)  £ 0)  4) 

£ 0£  M£  Q 

g P CP  0)  p H 

o 4 T3  I 
P CO  P (0  C 0)  >-l 
>P  -P  £ -P  3 O 84 


QQCAUZlO><(n84 

84 

84 

84 

(A 

84 

• 

• 

• 

• 

• 

• 

iH 

CN 

m 

fNfn^tinvor'OOO't-4 

rH 

rH 

12 


NSWC/WOL/TR  76-7 


Thus,  in  order  to  estimate  the  strength  degradation  of  composites  as 
a result  of  long  term  exposure  to  the  environment,  specifically  to 
humidity,  at  least  two  things  must  be  determined: 

1.  The  strength  (or  modulus)  retention  as  a function  of 
temperature  and  the  amount  of  moisture  absorbed. 

2.  The  rate  of  moisture  absorption  (diffusion)  as  a function 
of  temperature  and  relative  humidity.  Knowing  the  rate  should  allow 
one  to  predict  how  long  it  takes  to  reach  a certain  change  in 
strength. 

Of  course,  different  strength  parameters  such  as  tensile 
strength,  flexural  strength,  shear  strength,  compressive  strength 
and  their  corresponding  moduli  are  affected  differently  for  the  same 
change  in  moisture  concentration.  Also,  it  is  assumed  that  a 
certain  moisture  content  in  the  composite  corresponds  always  to  only 
one  strength  (modulus)  value  independent  of  the  way  the  moisture  was 
allowed  to  be  absorbed  in  the  seunple.  There  is  experimental 
evidence  that  the  strength  changes  with  moisture  pick-up  are  the 
same  whether  the  seunple  was  moisture  loaded  by  a water  boil  test  or 
by  water  soak  at  ambient  temperature  [5] . This  fact  might  be 
expected  as  long  as  the  moisture  concentration  distribution  is  the 
same  in  the  interior  of  the  composite.  The  same  overall  moisture 
concentration  might  also  correspond  to  different  interior  distri- 
butions, in  which  case  there  may  or  may  not  be  only  one  strength 
value  for  a given  value  of  percent  moisture. 

The  time  to  reach  a certain  moisture  concentration  (or  better, 
moisture  distribution)  is,  of  course,  dependent  on  the  composite 
geometry  and  the  temperature  and  humidity  profile,  but  can  be 
calculated  in  principle  if  the  diffusion  coefficients  (and  their 
temperature  and  concentration  dependences)  are  known. 

Thus,  one  of  the  very  important  parameters  to  be  measured 
for  the  prediction  of  property  changes  in  composites  is  their 
effective  diffusion  coefficients  of  moisture  along  the  principal 
axes  of  diffusion.  Of  special  interest  is  the  diffusion  coefficient 
across  the  fiber  direction  because  of  the  sheet  geometry  of  most 
composites. 

C.  Methods  For  Determining  Diffusion  Coefficients  in  Polymers 

The  two  most  widely  applicable  methods  for  measuring  diffusion 
coefficients  (D)  are  based  on  either  steady  state  flow  rate 


[5]  Personal  communication  with  N.  Judd  Royal  Aircraft  Establishment, 
Farnborough,  England,  and  Hercules,  Inc. 
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determinations  through  a membrane  or  on  the  analysis  of  sorption 
data  [6  to  3 ] . 

The  flux  (J)  through  a membrane  is  described  by  equations  la 
and  lb. 


J 


-D 


dc 
c dx 


(la) 


/ 


Jdx  = 


(lb) 


where  = the  diffusion  coefficient  (which  may  depend  on 
concentration),  c = concentration,  x = distance  through  the  membrane. 


In  a steady  state  condition  the  flux  is  independent  of  x and 
equation  (2)  follows. 


J 


D dc 
c 


(2) 


The  dependence  of  D on  concentration  can  be  determined  by 
measuring  the  flux  at  different  concentrations  on  one  side  of  the 
membrane  while  the  other  side  is  kept  essentially  at  zero 
concentration . 

The  other  method  for  determining  D is  based  on  absorption  or 
desorption  measurement  of  the  diffusant  in  plate,  sphere  or  rod 
samples . 

The  solution  of  Pick's  second  equation  of  diffusion  (3) 


Tt 


= D 


.2  .2  1 
3c.oC.oc 
— 7 — 7 — 7 

3x  3y  3z 


(3) 


[6]  J . Crank , "The  Mathematics  of  Diffusion",  Oxford  Univ.  Press. 
1956. 

[7]  J.  Crank  and  G.  S.  Parks,  "Diffusion  in  Polymers"  Acad.  Press 
(1968)  . 

[8]  G.  J.  Van  Amerongen,  "Diffusion  in  Elastomers,"  Rubber  Chem. 
and  Techn.  37,  1065,  (1964). 
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for  a two-dimensional  plate  with  infinite  extension  and  exposure  to 
constant  concentration  Ci  and  C2  of  the  diffusant  on  either  sides 
and  with  constant  initial  internal  concentration  Cq  can  be  expressed 
by  the  series  solution  (4)  [6]  . 


1 8_ 


m = o (2m  + 1) 


^ exp  [-D(2m  + l)^n^t/£^]  (4) 


where 

= Weight  increase  after  time  t,  = Weight  increase  at 
saturation  level,  = thickness  of  plate. 

For  Fickean  behavior  of  diffusion  the  plot  of  vs. 

is  practically  a straight  line  up  to  a value  of  0.6  of 
Thus,  the  initial  absorption  or  desorption  curves  of  uniformly 
loaded  plate  specimens  can  be  used  to  determine  D by  either 
expression  (5)  (p.  241  of  [6])  or  (6). 


D = 0.049/(t/i^) 


(5) 


where  the  subscript  *5  indicates  that  t is  taken  as  the  time  for 
to  reach  the  value  0.5. 


(6) 


If  the  initial  gradient,  G = d /d  ( t/H^ ) ^ ^ is  observed  in 

a sorption  experiment  in  which  D is_concentration-dependent , then 
the  average  diffusion  coefficient,  D,  deduced  from  (6)  is 


We  have  used  this  expression  to  calculate  the  diffusion 
coefficients  from  the  sorption  measurements. 

D.  Diffusion  of  Moisture  in  Fiber  Reinforced  Composites 

Since  composites  are  anisotropic  media  Fick's  second  law  of 
diffusion  takes  a more  complicated  form  (8)  than  for  isotropic 
media  (3).  The  diffusion  coefficient  has  a tensor  form. 
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3c  , n 3^C  ^ n 9^C  ^ r^  9^C  ^ ^ ^ 9^C 

at  T:T  °22  ^ ^33  ^ '^23  °32^  3y  &z 


3x 


^°31  ^13^  az'^Sx  ^°12  '*’  °21^  &x^&y 


(8) 


where  x,  y and  z are  Cartesian  coordinates. 

Equation  (8)  can  be  transformed  into  one  with  three  principal  axes 
of  diffusion,  with  principal  diffusion  coefficients  , 02/ 
and  D3  and  (8)  can  be  expressed  by  (9) 


3c 

3t 


(9) 


where  C/  n#  C/  are  the  principal  axes.  Various  symmetry  consi- 
derations such  as  found  in  crystallographic  systems  may  lead  to  still 
simpler  expressions. 


In  a system  with  more  than  one  solid  component,  such  as  a 
filled  plastic  or  a fiber  reinforced  composite,  the  overall  diffu- 
sion is  obviously  dependent  on  the  diffusion  in  both  substances,  on 
the  equilibrium  distribution  of  diffusant  in  substance  A and  B as 
expressed  by  the  distribution  coefficient  k = C^/Cb,  and  on  the 
spatial  arrangement  of  A and  B in  the  solid  phase. 

E.  Effective  Diffusion  Coefficient  in  Unidirectional  Carbon 
Fiber  Composites 

From  the  symmetry  of  unidirectional  carbon  fiber  composites 
two  of  the  three  principal  diffusion  coefficients  are  shown  to  be 
equal  (see  Figure  la) . The  question  to  be  answered  is,  how  do  the 
rates  of  diffusion  change  with  respect  to  the  unfilled  matrix  as  we 
change  the  fiber  content  and  the  fiber  direction. 

If  the  diffusion  of  moisture  in  the  carbon  fiber  is  negligible 
compared  with  that  in  the  resin,  the  problem  is  simplified.  A 
further  assumption  is  that  there  are  no  voids  in  the  resin  or 
between  resin  and  fiber  and  that  there  is  no  significant  mass 
transport  along  the  resin/fiber  interface. 

In  compacting  fibers  in  a resin  under  pressure  there  is  a 
higher  probability  for  the  fibers  to  assume  a hexagonal  packing 
(Figure  lb)  than  a tetragonal  arrangement  (Figure  Ic) , although  in 
a real  composite  a mixed  geometry  due  to  many  latice  defects 
probably  exists. 
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The  reduced  diffusivity  parallel  to  the  fiber  direction  in  a 
carbon  fiber  epoxy  composite  with  hexagonal  fiber  arrangement  is  the 
ratio  of  the  effective  diffusion  coefficient  to  the  diffusivity  of 
the  unobstructed  resin. 


e H 


h = 1 - 


V, 


(14) 


For  a tetragonal  fiber  arrangement  one  can  arrive  at  similar 
relations  (from  Figure  Ic)  as  indicated  below 


nr' 


(2r  + d) 


(15) 


2r'V  d = 


(16) 


M ^ 1.7724 

V/  t Vj 


(17) 


(18) 


The  fiber  volume  fraction  for  the  highest  packing  density 

(d=0)  for  composites  with  hexagonal  fiber  arrangement  is  ttV3/6 
and  with  tetragonal  arrangement  it  is  tt/4  . 

The  ratios  of  d/r  for  the  hexagonal  and  tetragonal  fiber 
arrangements  as  a function  of  percent  fiber  volume  is  shown  in 
Figure  2. 

Effective  diffusion  coefficients  for  diffusion  of  water  vapor 
across  the  fiber  direction  in  a composite  material  with  a hexagonal 
or  tetragonal  arrangement  of  fibers  may  be  determined  numerically 
by  solving  a Poisson  equation  (with  appropriate  boundary  conditions) 
on  a unit  cell  of  the  material  (a  detailed  demonstration  that  it 
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To  simplify  the  discussion,  we  consider  a two-dimensional 
plate  through  which  moisture  diffuses  under  steady  state  condi- 
tions. If  there  are  no  fibers  in  the  matrix  the  amount  of  material 
diffusing  through  the  plate  is  given  by  Pick's  first  law 

Jr  = -Dr  ^ (10) 

2 

where  J = flux  in  g/cm  sec  through  the  surface  vertical  to  the 
flow  difection,  D^.  = average  diffusion  coefficient  of  the  resin, 

AC  = difference  in  concentration  of  the  diffusant  on  both  matrix 
surfaces,  I = thickness  of  the  sheet. 

If  fibers  are  introduced  into  the  matrix  parallel  to  the  flow 
direction,  the  effective  cross  section  for  the  diffusion  is  reduced 
proportionally  to  the  cross  sectional  area  of  the  introduced  fibers, 
i.e.,  proportionally  to  the  fiber  volume  fraction,  V^. 

The  cross  sectional  area  of  the  fibers  with  respect  to  a unit 


cell  of  a hexagonal  fiber  arrangement  is 


/2  (where 


r = radius  of  fiber) . The  total  cross  sectional  area  of  the  unit 

2 

cell  is  V,  = (2r  + d)  (sin  60)/2.  The  fiber  volume  fraction  is 

therefore^ 


^f  = ilTTeo 


= 3.6276 


(2r  + d) 


(2r  + d) 


The  volume  fraction  of  the  resin  matrix  is  Vj,  = 1 - V^.  The 
subscripts  h or  t refer  to  hexagonal  or  tetragonal  fiber  arrange- 
ments, and  the  subscript  f and  r refer  to  fiber  and  resin 
respectively . 

From  (11)  follows 


. 52503 


1.9046 


where  d « shortest  distance  from  fiber  to  fiber  surface. 
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is  sufficient  to  solve  only  within  a unit  cell  is  given  in 
Appendix  A).  Unit  cells  for  computing  the  effective  diffusion 
coefficient  of  water  vapor  perpendicular  to  the  fibers  in  a hexa- 
gonal and  in  a tetragonal  array  are  shown  in  Figure  1. 

To  determine  the  effective  diffusion  coefficient  dg  (in  the  y 
direction)  with  either  a hexagonal  or  tetragonal  arrangement  of  the 
fibers,  we  consider  the  following.  Allow  no  diffusion  across  the 
vertical  sides  of  the  unit  cell.  Let  the  bottom  and  top  faces  of 
the  unit  cell  be  exposed  to  uniform  concentrations  of  water  vapor 
so  that  the  water  vapor  concentration  in  the  matrix  along  the 
bottom  face  is  a constant  C2  and  so  that  the  water  vapor  concen- 
tration in  the  matrix  along  the  top  face  is  a constant  . A 
steady  state  concentration  water  vapor,  C(x,y),  will  then  become 
established  in  the  unit  cell.  The  effective  diffusion  coefficient 
dg  is  given  by 


dg  = -q  • y/(C^  - (r  + d/2)  (19) 

where  y is  the  height  of  the  unit  cell,  and  q is  the  flux  through 
any  y=constant  surface  in  the  unit  cell  (q  being  independent  of  y 
in  the  steady  state) . Letting  d (x,y)  E (the  diffusion 
coefficient  of  the  fiber)  if  (x,y)  is  a point  lying  in  the  fiber, 
and  letting  d (x,y)  e D^.  (the  diffusion  coefficient  of  the  resin) 
if  (x,y)  is  a point  in  the  resin,  one  has  that  q is  given  by 


r + d/2 

q = f - 3"  (x,y)  Cy(x,y)dx  for  any  yc  [o,y]  . (20) 


The  specific  differential  equation  and  boundary  conditions 
satisfied  by  c(x,y)  will  now  be  given.  There  is  a constant  )<: 

(the  distribution  coefficient,  which  depends  on  the  particular 
materials  in  the  resin  and  in  the  fiber)  so  that  at  any  point  p 
lying  on  an  interface  separating  the  resin  and  fiber  regions,  one 
has 


Cf(p)/Cr(p)  = k, 

where  C£(p)  (C  (p) ) denotes  the  limit  at  p of  the  concentration 
as  one  approaches  from  within  the  fiber  (resin) . Let  the  function 
u(x,y)  be  given  by 


c(x,y),  for  points  (x,y)  in  the  resin 
c(x,y)/k  for  points  (x,y)  in  the  fiber 
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so  then  u(x,y)  will  be  continuous  across  any  interface  between 
matrix  and  fiber.  If  p is  a point  on  an  interface  between  resin 
and  fiber,  let  n denote  a unit  normal  to  the  interface  curve  at  p. 
Then  the  condition  of  continuity  of  flux  at  p is 


-D-  • k • u^(p)  = -D  u^  (p)  (22a) 

f n r n 

f r 

where  ^’^^^(p))  denotes  the  normal  derivative  of  u at  p calcu- 

lated from  within  the  fiber  (resin) . In  the  case  = 0 the  left 
side  of  (22a)  is  defined  to  be  0.  The  differential  equation  and 
other  boundary  conditions  satisfied  in  the  unit  cell  by  u are; 


0 = D Au  = D (u  + u ) in  the  resin,  (22b) 

r r XX  yy 

0 = • kAu  in  the  fiber  (22c) 

with  the  understanding  that  when  D-  = 0 conditions  on  u holding 
in  the  fiber  region  are  to  be  omitted, 

u(x,o)  = C2  and  u(x,y)  = for  0<x<r+d/2,  (22d) 

u (o,y)  = 0 and  u (r+d/2,y>  * 0 for  0 y y (22e) 

X X 


Note  that  (20)  in  terms  of  u is 
r+d/2 

q = y*  -d(x,y)  Uy  (x,y)dx  for  any  yc  [o,y] , (23) 

o 


where  d(x,y)  = in  the  resin  and  d(x,y)  = k • in  the 
fiber . 

Numerical  approximation  of  (19) , (22) , and  (23)  was  done  as 
follows  (a  completely  detailed  discussion  is  given  in  Appendix  1) . 
Rather  than  solve  the  interface  problem  (22)  directly,  which  would 
involve  dealing  directly  with  the  geometry  of  the  curved  interface 
between  matrix  and  fiber,  a variable  diffusion  coefficient 
problem  (24-26)  approximating  (22)  was  solved  on  the  rectangular 
unit  cell; 


u(x,o)  = C2 
Ux(o,y)  = 0 


and  u(x,y)  = Cj^ 
and  u^  (r+d/2, y) 


for  0<x<r+d/2 
= 0 for  0<y<y 
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V*d  (x,y)Vu  =0  on  0<x<r+d/2,  0<y<y  (26) 

where  away  from  the  interface,  d = D^.  in  the  resin  and  d = k*Df 
in  the  fiber,  while  in  a narrow  zone  near  the  interface  d changes 

fpntinuously  from  D^.  to  k-D^.  In  the  situation  Df  = 0 , 

was  defined  to  be  k/1000  in  the  fiber,  and  varied  from  D^.  to  k/1000 
near  the  interface. 

The  Equations  (24-26)  were  solved  using  the  D'Yakonov  ADI 
(alternating  direction  implicit)  finite  difference  method  on  a grid 
with  uniform  mesh  spacings  Ax  and  Ay  in  the  x and  y directions, 
respectively.  The  flux  q given  by  (23)  (with  d(x,y))  replaced  by 
d(x,y))  was  evaluated  by  differencing  u to  obtain  Uy  and  by  using 
the  trapezoidal  quadrature  rule  to  evaluate  the  integral. 

F.  Thermal  and  Electrical  Analogs  to  Diffusion  Processes 

A thermal  analogue  to  the  diffusion  of  moisture  through  a 
medium  with  cylindrical  obstacles  perpendicular  to  the  flow  direc- 
tion has  been  treated  by  G.  S.  Springer  and  S.  W.  Tsai  [9] . These 
authors  have  calculated  the  thermal  conductivity  of  a system  with 
tetragonal  arrangement  of  cylindrical  rods  embedded  in  a matrix 
with  different  thermal  conductivity.  They  obtained  the  following 
equation 


tan 


_j^Jl  - (B'^Vj/n)\ 
' 1 + ^B^V^/n 


B = 2 ^ - 1^  (27) 

^22  “ conductivity  perpendicular  to  the  axes  of  the 

cylinders,)^  = thermal  conductivity  of  the  matrix,  = volume 
fraction  or  the  cylindrical  rods. 


If  kj  becomes  zero  (or  for  the  analog  in  diffusion,  the  diffusion 
coefficient  of  the  fiber  becomes  zero) , the  second  term  goes  to 
zero. 


At  equilibrium  the  difference  in  concentration  of  diffusant 
in  phase  (f)  and  phase  (r)  can  be  expressed  by  the  ratio 


[9]  G.  S.  Springer  and  S.  W.  Tsai,  "Thermal  Conductivity  of 
Unidirectional  Materials"  J.  Comp.  Mat.  166  (1967) . 
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k = c^/Cj.  (the  distribution  coefficient  between  resin  and 
fiber).  If  Henry's  law  applied,  the  analog  kDf/Dj-  can  be 
substituted  for  kf/kr  so  that  the  Springer-Tsai  equation 
becomes 


Another  analog  has  been  calculated  for  an  electrical  conductivity 
problem  by  Lord  Rayleigh  [10]  and  Runge  [11]  . Again  the  system  was 
composed  of  tetragonally  arranged  cylinders  in  a matrix  with  different 
conductivity.  An  equation  was  derived  which  gave  the  overall  conduc- 
tivity in  terms  of  the  conductivities  of  the  cylinder  and  matrix 
materials 


2V, 


1 + 


1 - 


K, 


1 - 


m 


m 


K 


+ - 


m 


1 + 


■'m 


X 0.3058  - 


(29) 


which  reduces  to  equation  (30)  when  the  conductivity  of  the  cylinders 
become  zero 


1 - 


2V, 


1 + - .3058 


(30) 


The  conductivity  in  the  Rayleigh  equation  does  not  become  zero  as 
it  should  for  the  highest  packing  density  of  tetragonal  cylinder 
arrangement  (see  Figure  3) . No  explanation  was  given  by  the  author. 

There  is  a considerable  difference  at  the  lower  fiber  volume 
fractions  between  the  Rayleigh  curve  and  that  of  Springer-Tsai  for 
zero  thermal  conductivity  of  the  cylinders. 


[16] Lord  Rayleigh,  Philosoph.  May.  3^,  481  (1892). 
[11] I.  Runge,  Z.  Tech.  Phys.  6,  61  (1925). 
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The  curve  we  obtained  from  the  finite  difference  approach 
was  in  excellent  agreement  to  that  of  Rayleigh  (for  the  tetragonal 
fiber  arrangement)  up  to  about  .70  fiber  fraction.  Above  this 
density  the  Rayleigh  equation  fails.  The  reason  why  the  Springer- 
Tsai  model  gives  lower  values  of  reduced  dif fusivities  over  the 
whole  range  is,  that  it  is  based  on  a simple  thermal  model  where 
the  composite  elements  are  connected  in  series  [9] . This  means 
that  a single  fiber  running  through  the  matrix  (vertical  to  the 
flow  direction)  has  the  same  effect  as  a whole  sheet  of  graphite 
with  the  sam.e  thic)<ness  as  the  fiber  being  placed  in  the  flow 
direction . 

That  there  is  a difference  can  be  shown  by  the  same  finite 
difference  method  if  the  following  boundary  value  problems  are 
solved : 


c 

a) 

$ 

.V 

b) 

for  the  unobstructed 
region  dx(2r+d) 

for  the  region  within 
(r+d)  X (2r+d) 


If  this  simple  thermal  model  would  be  applied  for  the  hexagonal 
fiber  arrangement  it  would  fail  already  at  tt/8  \/3  fiber  volume 
fraction  because  this  model  does  not  talce  into  account  that  the 
diffusant  can  flow  around  the  fibers.  The  difference  of  various 
models  is  given  in  Table  3 and  Figure  3. 


G.  Diffusion  Coefficient  in  an  Arbitrary  Angle  to  the  Fiber 
Direction  of  a Unidirectional  Carbon  Fiber  Composite 

The  principal  axes  of  diffusion  are  parallel  and  perpendicular 
to  the  fiber  direction.  D^ , therefore,  constitutes  a rotational 
ellipsoid  with  two  principal  diffusion  coefficients  being  equal. 

For  the  general  case  one  can  state  that  the  diffusion  coefficient, 

D,  at  right  angles  to  surfaces  whose  normals  have  the  direction 
cosines  m,  n relative  to  the  principal  axes  of  diffusion  is  given 
by 


D 


l^D  + 
x 


2^ 
m D 


n^D 


(31) 
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Table  3:  Reduced  Dif f usivities  as  a Function  of  Fiber  Volume  Fraction 
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and  the  corresponding  fluxes  are 


J 


X 


-D 

X 


9c 

9y 


J 

z 


-D 

z 


9c 

9z 


Thus,  knowledge  of  Dg||  and  DgX  is  all  that  is  required  to 
describe  diffusion  in  a unidirectional  carbon  fiber  composite. 

For  instance,  the  reduced  diffusivity  in  a unidirectional 
carbon  fiber  composite,  which  is  rotationally  symmetric  in  the 
fiber  direction,  can  be  expressed  as  follows. 


D 


ea 


2 2 
D^lj  cos  a + sin  a 


(32) 


H . Effective  Diffusion  Coefficients  in  Cross-Plied  Laminar 
Composites 

If  a composite  is  made  of  several  layers  of  the  same  prepreg 
but  with  the  layers  laid  up  such  that  the  fiber  directions  are 
different  for  each  layer,  there  is  basically  no  difference  in  the 
effective  diffusion  coefficient  perpendicular  to  the  fiber  direction 
if  the  average  fiber  volume  fraction  is  known.  The  reason  is  that 
a laminar  ply  consists  usually  of  many  layers  of  fibers  in  the  same 
direction.  The  effective  diffusion  coefficients  in  various  direc- 
tions in  the  plane  of  the  sheet  is  a superposition  of  all  the 
single  plies  and  their  respective  fiber  directions.  Each  ply  can 
be  treated  as  a unidirectional  composite  for  which  Dg|  and 
Dg  are  known  as  a function  of  their  fiber  volume  fraction.  Of 
course,  the  angular  dependence  of  the  diffusion  rates  in  each 
ply  is  not  independent  of  the  other  plies  since  there  is  an 
additional  requirement.  In  order  to  maintain  equilibrium,  the 
concentrations  at  the  bounaries  between  plies  (m  and  n)  must  be 
continuous.  For  instance,  if  the  flow  is  parallel  to  ply 
m ( = 0®),  then  the  flow  in  ply  n 0®)  is  increased  because  there 
is  an  outflow  from  m into  n which  has  a slower  rate  of  diffusion 
in  this  particular  direction.  Thus,  for  the  same  reason,  the 
effective  diffusion  rate  in  ply  m (in  the  0®  direction)  is  retarded 
as  compared  to  the  original  rate  of  diffusion  parallel  to  the 
fibers  without  neighboring  plies  (with  different  fiber  directions) . 

I.  Measurements  of  the  Diffusion  Coefficients  (DgX^  Carbon  Fiber 
Reinforced  Composite 

The  effective  diffusion  coefficients  (Dg^)  of  6 carbon  fiber 
composite  systems  (made  of  3 resins:  Narmco  5208,  DER  332/DADS 
and  Epon  1031/NMA  were  determined  as  a function  of  temperature  and 
relative  humidity. 


25 


NSWC/WOL/TR  76-7 


a.  Equilibrium  Concentration  of  the  Composites  as  a function 
of  temperature  and  relative  humidity. 

Calculations  show  that  it  would  take  more  than  100  years  for 
the  composite  plate  of  1 cm  thickness  to  reach  equilibrium  at  room 
temperature;  therefore,  the  equilibrium  concentration  was  measured 
on  thin  composite  shavings.  The  results  thus  obtained  compared 
well  with  calculated  values  from  resin  equilibrium  measurements, 
also  obtained  from  shavings  or  films.  The  calculations  were  made 
with  the  assumption  that  only  the  resin  absorbs  moisture.  The 
constant  humidity  exposures  of  30%  and  80%RH  were  carried  out  by 
adjusting  the  humidity  with  various  saturated  salt  solutions 
(Table  4)  for  different  temperatures.  The  results  are  shown  in 
Table  5.  Figure  4 shows  the  equilibrium  concentration  profile  of 
Narmco  5208  as  a function  of  temperature  and  relative  humidity. 

The  equilibrium  concentration  in  the  samples  increase  with 
increasing  RH  as  is  expected  from  Henry's  law,  which  states  that 
the  saturation  point  is  proportional  to  the  vapor  pressure.  However, 
in  some  cases  the  equilibrium  concentration  seems  to  decrease 
slightly  with  increasing  temperature  at  constant  relative  humidity 
in  spite  of  the  fact  that  at  constant  RH  the  absolute  concentration 
of  moisture  in  the  surrounding  air  increases  with  increasing 
temperature  (see  Figure  5) . The  reason  for  the  lower  equilibrium 
concentration  at  higher  temperature  may  be  due  to  the  fact  that 
the  surface  adsorption  is  lower  at  higher  temperature.  Thus,  the 
equilibrium  concentration  of  the  interior  seems  to  be  governed 
by  the  surface  adsorption.  Water  submersion  may  show  different 
results . 


b.  Fractional  Weight  Absorption  of  Carbon  Fiber  Composites. 

2 

5 X 5 cm  composite  plate  samples,  sealed  around  the  edges  with 
aluminum  foil  to  prevent  diffusion  parallel  to  the  fiber  direction, 
were  exposed  at  various  constant  temperatures  and  constant  RH 
values.  The  fractional  weight  increases,  M^/Moo,  (where 
is  the  percent  weight  increase  after  the  exposure  time  t and  Moo 
is  the  percent  equilibrium  concentration  as  determined  under  [a] ) , 
were  plotted  vs.  t^  for  the  composite  samples  Epon  1031/HMS, 
at  30%  and  80%  RH  at  30®,  45®,  60®  and  75®C  (See  Figures  6 and  7). 

c.  Effective  Diffusion  Coefficients  (De_|_)  of  Carbon  Fiber 
Composites  as  a Function  of  RH  and  Temperature. 

From  the  slope  of  the  sorption  curves  the  diffusion 
coefficients,  Dg^  were  calculated  according  to  equation  7.  The 
results  are  given  in  Table  6 and  as  three-dimensional  projections 
in  Figures  8 through  13. 

Since  the  composites,  made  from  the  T300  and  HMS  fibers,  only 
differ  about  4 percent  in  fiber  volume  concentration,  the  difference 
in  diffusion  coefficients  between  the  two  sets  of  composites  is 
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Table  4:  Relative  Humidities  Above  Saturatgd  Salt 

Solutions  At  Various  Temperatures^ 


Temp. 

rc) 

Approx . 

RH  < 

'^>33 

< ^.55 

< 'v-80 

30 

MgCiL2 

(33) 

Na2Cr20^  • 2H2O 

(57) 

NH^Ci  (81) 

45 

MgCi.2 

(32) 

Na2Cr20^  • 2H2O 

(54) 

NH^C)l  (79) 

60 

MgCil2 

(30) 

Na2Cr20^  • 2H2O 

(55) 

KC£  (80) 

75 

MgCJl2 

(29) 

Na_CrO.  (56) 
2 4 

KC£  (80) 

(a)  The  numbers  in  parentheses  are  the  RH  values  reported  in  the 
literature  for  the  respective  temperatures. 
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Table  5 

Equilibrium  Concentration  of  Moisture  in  Epoxy 
Resins  at  Various  Temperatures  (at  30  and  80%  RH) 


Resin 


Exposure  Exposure  Equilibrium 

to  (S  RH)  at  (®C)  Concentration  (%) 


5208 

80 

30 

4.64 

5208 

80 

45 

4.57 

5208 

80 

60 

4.50 

5208 

80 

75 

4.50 

332 

80 

30 

2.84 

332 

80 

45 

2.61 

332 

80 

60 

2.48 

332 

80 

75 

2.07 

1031 

80 

30 

1.67 

1031 

80 

45 

1.44 

1031 

80 

60 

1.23 

1031 

80 

75 

0.90 

5208 

30 

30 

1.93 

5208 

30 

45 

1.83 

5208 

30 

60 

1.80 

5208 

30 

75 

1.75 

332 

30 

30 

1.31 

332 

30 

45 

1.27 

332 

30 

60 

1.04 

332 

30 

75 

0.82 

1031 

30 

30 

0.60 

1031 

30 

45 

0.58 

1031 

30 

60 

0.52 

1031 

30 

75 

0.44 

5208 

100 

30 

5.93 

5208 

100 

75 

6.10 

5208 

100 

100 

5.98 

The  equilibrium  concentrations  have  been  determined  on  resin 
shavings,  or  resin  films. 
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expected  to  be  smaller  than  the  experimental  error  of  the  measure- 
ment. However,  if  the  ratios  of  the  Dgj^'s  of  the  HMS/T300 
composites  (see  Table  7)  are  listed,  then  a consistently  higher 
Dg  is  observed  for  the  HMS  composites  of  the  same  resins  and 
the  same  exposure  conditions.  The  diffusion  coefficients  for 
Narmco  5208,  DER  332/DADS  and  Epon  1031/MNA  HMS  composites  are  on 
the  average  about  40%,  60%  and  20%,  respectively,  higher  than  the 
corresponding  T300  composites.  The  T300  fiber  had  an  epoxy  sizing 
for  better  handleability  while  the  HMS  fiber  had  none.  This 
treatment  of  the  T300  fibers  with  an  epoxy  surface  coating  appears 
to  lead  to  a better  wetting  of  the  fibers  by  the  resin,  therefore, 
to  a reduction  of  inter facial  void  formation.  Obviously,  the 
diffusion  rates  through  the  voids  would  be  very  much  larger  than 
if  the  same  volume  were  filled  with  resin.  The  measurement  of 
void  content  by  the  resin  digestion  method,  however,  is  not  sensi- 
tive enough  (if  differences  are  due  to  voids) . 

Very  little  is  Icnown  about  the  interface  in  carbon  fiber 
composites.  The  questions  of  whether  there  is  a significant 
diffusion  along  the  interface  has,  to  our  knowledge,  not  been 
determined.  This  interface  may  play  a role  similar  to  the  grain 
boundaries  in  metals  and  ceramics.  Low  temperature  diffusion  of 
He,  H2  or  argon  might  give  some  insight  as  to  the  microporosity, 
particularly  along  the  resin-fiber  interface. 

d.  Distribution  of  Moisture  in  Composites  as  a Function  of 
Thickness,  Time,  Relative  Humidity  and  Temperature. 

The  solution  of  Pick's  second  diffusion  equation  for  a sheet 
geometry  is  shown  in  Figure  14.  The  solution  is  given  in  dimension- 
less parameters:  Dt/h^,  the  coordinates  are  given  as 
C“Co/Cj^-C„  and  X/h,  where  D = diffusion  constant,  t = time 
h = is  half  the  thickness  of  sheet,  Cq  = initial  uniform  concen- 
tration, = concentration  kept  constant  at  the  surface.  This 
graph  applies,  therefore,  for  all  substances  following  Fickean 
diffusion  behavior. 

Figure  15  allows  very  quickly,  without  the  use  of  a computer, 
the  determination  of  the  moisture  distribution  inside  a solid 
material,  once  the  diffusion  coefficient  of  moisture  has  been 
determined. 

2 

After  integrating  the  Dt/h  curves  of  Figure  15,  a plot  of 
vs.  lDt/h2]l/2  (again  of  dimensionless  parameters)  can 
be  derived.  This  plot,  shown  in  Figure  15,  permits  the  deter- 
mination of  the  internal  distribution  of  the  diffusant  by  simply 
measuring  the  sample's  fractional  weight  increase,  M^/M„» 
assuming  D has  been  previously  determined.  The  measured  M^/Moo» 
corresponds  to  a certain  [Dt/h^]^/2  in  Figure  15  from  which 
Dt/h2  is  determined  and  interpolated  on  Figure  14  to  give  the 
internal  distribution. 
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Table  7 

A comparison  of  the  effective  moisture  diffusion  coefficients  (D  ) 
in  the  HMS  and  T300  carbon  fiber  composites. 


Exposure  Temp. 

»c 

RH 

% 

Ratio  of  D . (HMS)/D 

eL 

Narmco  5208  DER  332/DADS 

. (T300) 

el- 

Epon  1031/NMA 

30 

30 

1.88 

1.10 

1.05 

45 

30 

1.60 

1.51 

1.36 

60 

30 

1.56 

2.08 

1.29 

75 

30 

.98 

1.61 

.92 

30 

80 

1.50 

1.42 

1.11 

45 

80 

1.12 

1.48 

1.36 

60 

80 

1.27 

2.16 

1.12 

75 

80 

1.16 

1.50 

1.51 

Average 
Coeff.  Var 
Stand.  Deviation 
of  the  Mean 


1.38  1.61  1.21 

7.70  7.69  5.72 

.1063  .1238  .0693 
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Of  course,  a finite  difference  method  for  calculating  the 
distributions  by  means  of  a computer  can  always  be  employed.  There 
is,  for  instance,  the  Crank-Nicolson  method  described  in 
reference  [6] . 

For  illustrative  purposes  a Narmco  5208/T300  composite  is 
considered.  The  samples  of  1 cm  thickness  are  exposed  to 
a)  80%  RH  at  25®C,  b)  30%  RH  and  25®C,  c)  80%  RH  and  35‘>C, 
and  also,  d)  a plate  of  0.25  cm  thickness  is  exposed  to  35 ®C  and 
80%  RH. 

The  moisture  equilibrium  concentrations  at  25®C  are  1.26%  for 
80%  RH  and  .582%  for  30%  RH,  and  at  35®C  and  80%  RH  the  equilibrium 
concentration  is  1.25%.  The  diffusion  coefficients  (see  Figure  8) 
at  25“C  and  30%  and  80%  RH  are  1 x 10“^®  cm^/sec  and  6 x 10“^^ 
cm^/sec  respectively  and  at  35®C  and  80%  RH  the  diffusion 
coefficient  is  «*2.5  x 10“^^  cm^/sec. 

Using  Figure  14  the  isochronal  distribution  of  moisture  inside 
the  composite  can  be  obtained  for  various  times;  these  distributions 
are  shown  in  Figure  16.  From  these  graphs  it  is  apparent  that  the 
time  required  for  98%  saturation  differs  widely.  Saturation  times 
depend  on  the  square  of  the  thickness  of  the  composite  and  on  the 
temperature,  which  roughly  doubles  the  rate  for  every  10 ®C  increase 
in  temperature. 

The  percent  weight  increase  with  time  for  these  4 conditions 
can  be  easily  obtained  from  Figure  15  by  using  the  relation  of 
Mt/Mco  vs.  [Dt/h2]l/2^  The  results  are  shown  in  Figure  17. 

For  varying  boundary  conditions  with  changing  temperatures  and 
humidities  the  problem  is  more  complicated  but  can  be  solved  with 
a computer.  Details  of  how  to  predict  the  effect  of  real  outdoor 
weather  conditions  with  constantly  varying  boundary  conditions  will 
be  the  subject  of  a forthcoming  NSWC/TR. 

Two  more  distributions  for  the  DER  332/DADS-HMS  and  Epon  1031/ 
NMA-HMS  composites  are  shown  in  Figures  18  and  19.  Both  moisture 
distributions  were  calculated  for  a 1 cm  thick  plate  exposed  at 
25®C  at  80%  RH. 

e.  Temperature  Coefficients  of  Diffusion 

The  temperature  coefficient  of  moisture  diffusion  for  the 
three  HMS  fiber  composites  at  30%  RH  was  obtained  from  Figure  20 
where  the  logarithms  of  the  diffusion  coefficients  were  plotted  vs. 
1/T  (®K)  . Since  the  equilibrium  concentration  is  only  slightly 
dependent  on  temperature  (at  least  below  60°C)  , the  slopes  of 
these  straight  lines  might  be  considered  as  overall  activation 
energies  of  diffusion  of  moisture  across  the  fiber  direction.  The 
values  range  from  about  8 to  10  kcal/mole.  Although  the  slopes 
are  about  the  same,  there  is  a difference  of  an  order  of  magnitude 
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between  the  Epon  1031  and  the  Narmco  5208  composites.  The  diffusion 
rate  of  the  Epon  1031  composite  is  about  10  times  higher  than  the 
Narmco  5208  resin,  in  spite  of  the  fact  that  the  moisture  equili- 
brium concentration  in  the  1031  composite  is  only  1/3  of  the 
equilibrium  concentration  of  Narmco  5208.  Therefore  equilibrium 
moisture  absorption  is  not  necessarily  an  indication  of  the  rate  of 
diffusion.  Thus  one  ought  to  be  careful  in  judging  resins  simply  by 
their  equilibrium  moisture  absorption. 

Furthermore,  as  will  be  reported  in  a subsequent  technical 
report  on  the  mechanical  property  changes  on  these  composites,  a 
low  moisture  equilibrium  concentration  does  not  necessarily  mean 
that  the  mechanical  properties  changes  will  be  proportionally  small. 

f.  Sample  Preparation. 

Since  the  sample's  final  equilibrium  concentration  of  moisture 
(Moo)  depends  on  the  relative  humidity  of  the  surrounding  environ- 
ment, the  relative  humidity  at  which  the  samples  have  been 
prepared  makes  a difference.  At  the  surface  of  the  sheet  of  resin 
(or  composite)  the  saturation  level  is  reached  within  a very  short 
time.  This  is  governed  by  the  surface  adsorption.  This  may  be  the 
reason  why  the  rate  of  absorption  at  a given  temperature  is  different 
for  different  relative  humidities  (see  Figure  21).  Thus,  samples 
with  the  same  percentage  of  moisture  content  may  have  different 
distributions.  For  instance  from  Figure  21  one  can  see  that  (for 
a sample  of  0.370  cm  thickness  of  Narmco  5208,  at  30 ®C,  a 0.4% 
moisture  pick-up  is  reached  in  water  after  4.8  days,  at  80%  RH 
after  6.6  days,  and  at  30%  RH  after  15.5  days.  Obviously,  since 
the  boundary  conditions  are  different,  so  is  the  internal  moisture 
distribution. 

It  is  not  yet  known  how  the  mechanical  properties  of  composites 
vary,  if  at  all,  when  they  are  prepared  in  different  RH's  but  having 
the  same  total  moisture  concentration.  Such  differences  in  prepa- 
ration will  have  to  be  considered  for  sample  testing  designed  to 
predict  the  property  changes  in  various  climatic  environments. 

The  ultimate  equilibrium  concentration  of  moisture  inside 
the  composite  will  correspond  to  the  longtime  average  humidity  in 
the  environment.  The  perturbation  due  to  short-time  (daily)  humi- 
dity fluctuations  is_restricted  only  to  a thin  surface  layer  on 
the  order  of  x = 2 (HD  t )*s,  where  x is  the  depth  in  change 
of  moisture  concentration  due  to  the  environmental  humidity 
fluctuation,  D is  the  average  effective  diffusion  coefficient  and 
tQ  is  the  length  of  the  cycle.  For  Narmco  5208/T300  composite 
this  distance  would  be  0.005  cm  at  ambient  temperature  and  a 
humidity  fluctuation  around  80%  RH  over  a 24-hour  cycle. 

J.  Comparison  With  Experimental  Data 

There  are  not  many  numerical  data  published  in  the  literature 
that  would  permit  the  verification  of  the  calculated  data  with 
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experiments.  Hertz  [1]  has  measured  the  diffusion  coefficients 
Dg  II  and  Dei.  in  a carbon  fiber  reinforced  epoxy  composite 
(V f = 58.8%,  and  2.3%  voids).  The  measured  data  were  DqH 
= 6.34  X 10-8  cm^/sec  and  = 1.56  x 10"®  cm^/sec. 

Unfortunately  he  did  not  give  the  diffusion  coefficient  of  the 
neat  resin.  The  ratio  of  Dgn  /D^x  should  be  between  1.6  - 
1.7  (see  Table  3).  The  experimental  ratio  is  4.06,  indicating 
Dgii  is  larger  than  should  be  probably  because  of  elongated  voids 
in  the  fiber  direction  (perhaps  debonding  along  the  fibers) . 

Also,  the  measured  values  of  Dell  and  Dq^  are  surpris- 
ingly high.  From  Dg^  = 1.56  x 10~8  cm^/sec,  according  to 
Table  3 the  resin  should  have  a diffusion  coefficient  of  about 
6.5  X 10“8  cm^/sec,  which  is  rather  high  for  epoxy  resins. 

Thus,  even  low  void  contents  seem  to  have  a pronounced  effect  on 
the  diffusion  coefficient.  On  the  other  hand  the  diffusion 
coefficients  we  measured  on  the  Narmco  5208  composites  (see 
Table  6)  are  all  lower  than  the  values  calculated  from  the  finite 
difference  method.  The  HMS  fiber  composites  are  17-20%  lower  and 
the  T300  composites  between  36-53%  lower  than  the  values  calcu- 
lated for  a tetragonal  fiber  arrangement.  One  reason  could  be  that 
the  resin  samples  were  cast  from  a different  resin  batch  than  the 
resin  used  for  the  composite  preparation.  We  have  previously 
observed  on  a different  epoxy  resin  that  within  the  same  resin 
batch,  a difference  in  thermal  history  (for  instance  different 
gelation  temperatures  but  same  cure  and  post  cure)  give  a difference 
of  22%  in  resin  diffusion  coefficients.  Another  reason  could  be 
that  the  composite  equilibrium  values  was  determined  on  powder 
samples.  Since  these  fine  powders  can  equilibriate  very  rapidly 
there  is  the  possibility  of  a weighing  error  due  to  the  time 
required  for  removing  the  samples  from  the  exposure  chambers  and 
weighing  them. 


Conclusions  and  Recommendations 

A finite  difference  method  was  used  to  estimate  the  reduced 
dif f usivities  in  unidirectional  carbon  fiber  composites  with  hexa- 
gonal and  tetragonal  fiber  arrangement.  For  the  tetragonal  arrange- 
ment of  cylindrical  rods  these  derivations  were  compared  with  two 
models,  derived  by  other  investigators  for  the  analogous  electrical 
and  thermal  conductivities.  While  the  Rayleigh  model  gives  almost 
identical  results  with  ours  at  fiber  concentrations  below 
60  percent,  the  Springer-Tsai  model  deviates  considerably  over  the 
whole  range  of  fiber  volume  fractions. 

The  effective  diffusion  coefficients  of  6 carbon  fiber/epoxy 
composites  were  determined  (perpendicular  to  the  fiber  direction) 
as  a function  of  temperature  and  relative  humidity.  There  is  a 
reasonable  agreement  with  the  calculated  reduced  diffusivity  for 
the  .7  fiber  volume  fraction  of  these  composites. 
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The  effective  diffusion  constants  for  the  composites  follow 
the  order  0(332)^  ^(5208^  while  the  saturation  values 

follow  the  inverse  or3er  C,..  (5208)  (332)  > C,,  (io31).  This 

reversal  of  order  indicates  that  measurement  of  moisture  picltup 
alone  is  not  sufficient  for  determining  the  effect  of  moisture  on 
resins  or  composites. 

Diffusion  processes  should  be  further  investigated  since 
diffusion  measurements  seem  to  be  quite  sensitive  to  composite 
changes  (that  can  generate  micro  cracks)  and  which  might  not  be 
r^^adily  detected  by  mechanical  testing.  Furthermore,  diffusion 
measurements  should  be  combined  with  mechanical  and  thermal  cycling 
experiments . 

In  addition  to  moisture  di f f usivi ties , we  expect  that  low 
temperature  He,  Ar , or  II2  dif  f usivities  may  give  gome  additional 
information  about  the  resin-fiber  interface. 

A mathematical  model  should  also  be  derived  to  include  surface 
coatings  on  composites  to  determine  any  potential  benefit  of  these 
coatings . 
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FLOW(  J||i  PARALLEL  TO 
FIBER  DIRECTION 


b.  (HEXAGONAL  FIBER 
ARRANGEMENT) 


c.  (TETRAGONAL  FIBER 
ARRANGEMENT) 


FIG,  1 GEOMETRY  OF  UNIDIRECTIONAL  FIBER  REINFORCED  COMPOSITES. 
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TETRAGONAL  FIBER  ARRANGEMENT). 
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FIG.  3 A COMPARISON  OF  THE  REDUCED  DIFFUSIVITY  VERTICAL  TO  THE  FIBER 
DIRECTION  WITH  THERMAL  AND  ELECTRICAL  ANALOGS. 
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FIG.  4 EQUILIBRIUM  CONCENTRATION  OF  MOISTURE  IN  NARMCO  5208  RESIN  AS  A FUNCTION  OF 
TEMPERATURE  AND  RELATIVE  HUMIDITY 
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FIG.  5 CONCENTRATION  OF  MOISTURE  IN  AIR 
AS  FUNCTION  OF  TEMPERATURE 
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TIME  ’^2(h0URS) 

FIG.  7 FRACTIONAL  WATER  ABSORPTION  AT  VARIOUS  TEMPERATURES 
OF  EPON  1031/HMS  COMPOSITE  AT  80%  RH 
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fig.  8 OIFF.  COEFFICIENT  OF  WATER  FOR  NARMC0  5208/T300C-FIBER 
COMPOSITE  AS  A FUNCTION  OF  TEMP.  AND  REL.  HUMIDITY 
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FIG.  9 DIFF.  COEFFICIENT  OF  WATER  FOR  DER  332/T300  C-FIBER 
COMPOSITE  AS  A FUNCTION  OF  TEMP.  AND  REL.  HUMIDITY 
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FIG.  10  DIFF.  COEFFICIENT  OF  WATER  FOR  EPON  1031-T300C-FIBER 
COMPOSITE  AS  A FUNCTION  OF  TEMP.  AND  REL.  HUMIDITY 
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FIG.  11  DIFF.  COEFFICIENT  OF  WATER  FOR  NARMCO  5208/HMS  C FIBER 
COMPOSITE  AS  A FUNCTION  OF  TEMP.  AND  REL.  HUMIDITY 
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FIG.  12  DIFF.  COEFFICIENT  OF  WATER 
COMPOSITE  AS  A FUNCTION  01 
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TEMPERATURE  AND  RELATIVE  HUMIDITY 
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FIG.  14  CONCENTRATION  DISTRIBUTION  (MASTER  PLOT)  AT  VARIOUS  TIMES 
IN  THE  SHEET  -h<  X < h WITH  INITIAL  UNIFORM  CONCENTRATION  C„ 
AND  SURFACE  CONCENTRATION  C,.  NUMBERS  ON  CURVES  ARE  VALVES 
OF  DtV  . 


16  DISTRIBUTION  OF  MOISTURE  IN  A UNIDIRECTIONAL  5208/T300  EPOXY  COMPOSITE: 
a.  AT  25°C  AND  80%  RH  (THICKNESS  1 CM);  b.  AT  25°C  AND  30%  RH  (THICKNESS  1CM); 
c.  AT  35°C  AND  80%  RH;  d.  AT  36°C  AND  80%  RH  (THICKNESS  0.25  CM) 
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FIG.  21  RATE  OF  MOISTURE  PICK  UP  IN  A 0.370  CM  THICK  SHEET  OF  NARMCO  5208  RESIN  AT  30°  C 
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APPENDIX  A 

ALGORITHM  FOR  THE  NUMERICAL  CALCULATION  OF 
EFFECTIVE  DIFFUSION  COEFFICIENTS 

In  this  appendix  a complete  description  will  be  given  for  the 
numerical  method  used  to  calculate  effective  diffusion  coefficients 
perpendicular  to  the  fiber  direction  in  a carbon  fiber  composite. 

As  described  in  Section  E.  of  Discussion  and  Results,  the  geometry 
under  consideration  is  either  hexagonal  or  tetragonal  arrangement 
of  the  fibers  (Figure  1) . 

We  first  demonstrate  that  it  suffices  to  solve  the  diffusion 
problem  within  a unit  cell,  and  then  we  give  a detailed  description 
of  the  finite  difference  method  used  in  solving  (24-26) . 

Al.  Sufficiency  of  Solving  within  the  Unit  Cell 

It  suffices  to  show  that  the  solution  on  the  full  infinite 
plate  region  of  the  problem  can  be  generated  from  the  solution  on 
the  unit  cell,  and  that  (as  will  be  obvious  from  the  construction) 
the  effective  diffusion  coefficient  calculated  from  the  unit  cell 
is  identical  with  that  obtained  from  the  entire  domain  of  the 
problem . 

The  "wea)^"  formulation  of  (22)  is  (G  unit  cell)  ; 

/ (d(x,y)u  V + d(x,y)u  v )dxdy  = 0 (33) 

Jfc  ^ ^ y y 

for  all  square  integrable  test  functions  v having  square  integrable 
first  derivative  in  G (such  functions  v are  said  to  be  in  H' (G) ) and 
which  vanish  on  the  line  segments  y=0  and  y=y  for  0<x<  r + d/2.  The 
solution  u of  (33)  is  to  lie  in  H’ (G)  and  satisfy  (?25) . Because 
the  interface  is  smooth  and  meets  the  boundary  of  G at  right  angles, 
the  (unique)  solution  u of  (33)  has  square  integrable  second  deri- 
vatives within  the  resin  and  within  the  fiber  and  is  continuous  on 
G and  its  boundary  [12],  [13].  By  using  Green's  formula  on  (33) 

to  shift  derivatives  from  v onto  u,  one  has  (b"r+d/2, 

[12]  R.  bT  Kellogg , "On  the  Poisson  Equation  with  Intersecting 
Interfaces,  "Applied  Analysis,  4,  101,  (1975). 

[13]  D.  P.  Squire,  "On  the  Existence  and  Uniqueness  of 
Solutions  of  the  Poisson  Interface  Problem,  "Am.  J.  Math. 

85,  241,  (1963). 
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I E interface  curve,  n = unit  normal  to  I pointing  into  the 
resin)  ; 


x,y) Au- V 


+ 


(d(0,y)u^(0,y)v(0,y) 


d (b,y)Uj^  (b,y)  v(b,y) ) " ^r%^  ) 


(34) 


and  by  varying  v over  H' (G)  (with  v=0  when  y=0  and  when  y=y)  one 
recovers  (22) . 

Using  the  weak  formulation,  one  verifies  that  the  solution 
may  be  extended  horizontally  (i.e.  parallel  to  the  x-axis)  by 
reflection  (as  is  intuitively  obvious  from  the  symmetry  conditions 
(22e)).  To  extend  vertically,  we  use  the  Schwarz  reflection 
principle  [14] . By  addition  of  a constant,  we  may  take  C2=0 
in  (22d) , and  define 


u(x,  -y)  = -u(x,y)  for  0^y<y,  all  x 


(35) 


Similarly  one  may  continue  the  solution  an  arbitrary  number  of  cell 
heights  along  the  y-axis. 

A2 . The  D'Yakonov  ADI  Method 

Equations  (24-26)  have  the  form 

V*AVu  = (A(x,y)u  ) + (A(x,y)u  ) = 0 (36a) 

XX  y y 

on  the  rectangle  R = |o<x<b,  0<y<y| 

u(x,o)=C2  and  u(x,y)  = for  0<x<b  (36b) 

Ux(o,y)  = u^  (b,y)  = 0 for  0<y<y  (36c) 


[l4]  P.  R.  Garabedian,  "Partial  Differential  Equations, 
John  Wiley  and  Sons,  Inc.  1964. 
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where  for  (24-26)  A=d  and  b=r+d/2.  Following  [15],  we  now  describe 
the  D'Yakonov  ADI  (alternating  direction  implicit)  finite  difference 
method  as  it  applies  to  (36)  (i.e.  (24-26)). 

We  consider  (36)  to  be  the  steady  state  of  the  time  dependent 
problem 


= V»AVu  in  R 

u(x,o,t)  = and  u(x,y,t)  = 
u^(0,y,t)  = (b,y,t)  = 0 for 

u(x,y,t)  = u^(x,y)  in  R 


(37a) 

for  0<x<b,  t^  0 

(37b) 

0<y<y,  t^O 

(37c) 

(37d) 

where  Uq  is  given  initial  data.  By  choosing  Uq  and  solving  (37) 
one  may  approach  the  steady  state  solution  u(x,y)  of  (36). 

To  solve  (37),  a time  step  At  and  positive  integers  J and  K 
are  chosen;  the  spatia^  mesh  lengths  are  then  defined  to  ^e 
Ax  = b/(J-l)  and  Ay  = y/(K-l) . An  approximate  solution  U.  . of  37 

n 3 » 

is  then  obtained  (n=l,2,...)  where  U.  . refers  to  the  value  of  U 

J f ^ 

at  time  n*At  at  the  spatial  mesh  point  ( ( j-1) Ax , (k-1) Ay ) . By 
definition  U.  . = u„ ( ( j-1) Ax , (k-1) Ay)  . Let  r = At/  Ax^ , 

r^  = At/Ay^  and  use  the  notation  (for  any  variable  H defined 

at  the  mesh  points) 

“^x^  '^x  ^*j,k  ■ ^j  + l/2,k  (^j+l,k  ^j,k)  " ^j-l/2,k  ^^j,k  ^^j-l,k^ 

^ '^y  ,k+l/2  ( “ j ,k+l  ^j,k)“  ^j,k-l/2  ^^j,k  ^j,k-l^ 

where  k ~ * ((j+1/2)  Ax,kAy)  etc. 

The  D'Yakonov  formula  (i.e.  system  of  linear  equations)  for 
obtaining  from  U*?  . which  involves  intermediate  values  is; 


1 

2 


r 6 A6 

X X 


J 


vn+1 

j 


z"  , for  l<j<J  2<k<K-l 

3 » ^ — I 


(38a) 


[15]  A.  R.  Mitchell,  "Computational  Methods  in  Partial  Differential 
Equations,"  John  Wiley  Sons  Ltd.,  1969. 
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r 6 A6 

y y y 


) ^j!k  = 2<k<k-i 


(38b) 


where 


By  making  use  of  the  boundary  conditions  (37b,c)  we  define 


u"  , = = C,  for  0£j£J+l,  n=0,l,2,  ... 

J / X Z J 9^  X 


(40a) 


^o,k  ~ ^2,k,  ^J+l,k  ■ ^J-l,k  llk<K,  n 0,1,2, 


(40b) 


and  (also  by  symmetry) 


A 1 ~ ^ -i  ^7  j.  1 

2'  2'  ^ 


= A 


J - 1 


for  l<k<K 


,k 


(40c) 


and  from  (38b) 


''o,k  '^2,k,  J+l,k 


= Vj_j^  for  2£k^K-l  ,n=l , 2 , 


(40d) 


which  completes  the  Equations  (38)  . The  point  of  the  ADI  method 
is  that  (38a)  involves  K-2  separate  tridiagonal  linear  systems  of 
J equations,  rather  than  a total  system  of  J(K-2)  linear  equations 
(and  similarly  with  (38b))  which  results  in  a drastic  decrease  in 
the  work  necessary  to  solve  the  equations. 

A3 . Comments  on  Implementation 

Concerning  actual  computer  implementation  of  the  D'Yakonov 
method,  the  following  observations  are  useful.  Since  the  coeffi- 
cient A and  the  boundary  conditions  in  (37)  are  independent  of 
time,  the  decomposition  of  the  K-2  tridiagonal  linear  systems 
from  (38a)  (and  the  J systems  from  (38b))  can  be  done  once  and 
then  stored.  Letting 

Q = (1  + i r 6^  A 6^)  U°  for  0<j<J+l,  2<k<K-l 

jfiv  ^yy  y J f ^ — — 
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so 


Z°  , = (1  + ^ r 6 A6  )Q.  , for  l<j<J,  2<k<k-l 
],k  2 XX  X ],k  — — — 


one  notes  by  using  the  boundary  conditions  (37b,c)  which  gave  (40) 
that 


Q 


o , k 


^2,k,  °J+l,k 


, . for  2<k<K-l. 

w X ^ K 


(41) 


It  is  also  very  useful  to  note  the  following  method  for 
calculating  (n=l,2,...)  (16],  [17].  Suppose  we  have  just 

calculated  . . We  then  have  available  (assuming  proper  strategy 

in  the  computer  code).  , V.  . , and  z"  . . We  observe  that 

J f K J f ^ 3 f ^ 

from  (38b)  (and  dropping  the  subscripts  jrk) 
u"  - v"  = J ry  6y  A6y  u" 

from  which 


(1  + j r^6y  A5y)u'^  = 2U”  - v" 


so 


z"  = (1  + 1 '^x'^x^'^x^ 


but  by  (38a) 


from  which 


1 

2 


r 6 
X X 


AS  v'^ 

X 


[16]  S.  Leventhal,  personal  communication 

(17]  M.  Ciment  and  S.  Leventhal,  "Higher  Order  Compact 
Implicit  Schemes  for  the  Wave  Equation,"  Math.  Comp.,  29, 
985,  (1975). 
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and  so 


Z = 2U  - V + -r  6 A6  2U  + Z - V , and 
2 X X X 

= 2U^  - 2V^  + + r 6 A6 

XX  X 


(42) 


which  greatly  simplifies  the  computation  of  Z^  at  each  step. 
A4 . Calculation  of  Diffusion  Coefficients 


The  D’Yakonov  ADI  scheme  was  used  to  solve  (24-26)  by  relaxing 
the  corresponding  time  dependent  problem  toward  the  steady  state. 
Then  dg  was  obtained  using  (19)  and 


r+d/2 


q = 


/ 


d(x,y)  (x,y)dx 


for  any  ye  [o,y] 


(43) 


Since  actually  it  is  the  ratio  d /D  that  is  of  interest,  we  may 

e r 

scale  distance  so  that  b=l  and  time  so  that  D^=l.  For  convenience 

we  also  set  Cj^  = 0 and  C2=y  so  that  (from  (19))  d /D  = g. 

e r 

The  explicit  definition  of  d(x,y)  for  the  tetragonal  unit  cell 
is  as  follows.  A parameter  z>0  is  given  (in  general  we  took 
z = .02  with  J=K=21) . Given  (x,y) , let  p = ((x-b)2  + y2jl/2 

and  define  D^  = max(.0^1,  Df ) . _If  p r + z,  then  d(x,y)  = D ; 
while  if  r-z,  then  d(x,y)  = Dj*  k.  For  r-z^p^r+z, 


d(x,y)  = D^*k  + (p-(r-z))(D^  - k)/(2z). 


For  the  hexagonal  unit  cell  one  also  similarly  deals  with  the  fiber 
centered  at  (0,y).  In  the  hexagonal  case  we  generally  choose 
z=.02  with  J=21  and  K=41. 


The  calculation  of  (43)  was  done  at  several  different  y, 
values  (y^  ; (m  - 1/2) Ay)  using  the  trapezoidal  rule; 


m 


q^q 


n 

- = E 

i=l 


d(  (i-1)  Ax,y  ) (U*!*  , - , ) 

m 1-1,  m+1  1-1, m 
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Comparing  for  several  values  of  m j)rovidcd  a good  indji^cation  of 
when  n was"'suf f iciently  large;  once  uniformity  of  the  values 

was  obtained  (for  several  m) , further  increases  in  n produced 
little  change  in  the  flux  values. 

The  choice  for  initial  values  was  Uq(x,'/)  = y-y.  Other  choices 
of  Uq  attempting  to  better  approximate  u(x,y)  produced  some 

instances  of  erratic  and  slow  convergence  of  U.  P toward  the  steady 

J / K 

state.  The  program  was  written  so  that  a sequence  of  several 
different  time  steps  could  be  run  during  the  course  of  a single 
problem  in  order  to  obtain  faster  approach  toward  steady  state 
values.  With  Df  = 0,  k=2,  tetragonal  geometry,  J=K=21,  z=.02  and 
fiber  volume  fraction  in  the  range  0 - .7,  an  appropriate  sequence 
of  time  steps  was  found  to  be  40  steps  with  At  = .05,  followed 
by  20  steps  with  At=20.0,  followed  by  40  steps  with  At=.5. 

With  Dp=0,  k=2,  hexagonal  geometry,  J=21,  K=41,  z=.02,  and  fiber 
volume  fraction  in  the  range  0-.85,  an  appropriate  sequence  of 
time  steps  was  found  to  be  60  steps  with  At=.02,  followed  by 
60  steps  with  At=20.0.  Convergence  tended  to  be  more  rapid 
with  positive  values  of  Df  in  (0,1).  These  choices  of  J,K,  At,  z 
were  found  by  observing  the  effects  of  varying  these  parameters  on 
the  values  of  dg  thus  obtained. 

A5 . A Numerical  Example 

We  give  results  obtained  from  Df  = 0 , k=2,  hexagonal  fiber 
geometry,  J=21,  K=41,  z=,02,  fiber  volume  fraction  =.4,  and  the 
choices  of  At  stated  above.  Fluxes  were  calculated  across 
y„,  = (m- 1/2)  Ay  for  m=l,  4,  10,  20,  40.  The  corresponding  flux 
values  at  t=o  (i.e.  for  Uq)  were  .3341,  .3554,  .4802,  1.0,  .3341. 
After  60  steps  with  At=.02  the  flux  values  were  .4279,  .4280, 

.4292,  .4321,  .4279.  After  60  more  steps  with  At=20.0  all  five 
flux  values  were  .4289. 
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Glossary 


Symbols 


C , c 


d 


D,d 


D 


r 


D 


f 


D 


D 


e' 


d 

e 


h 


m 


m 


q 


M 

00 


r 

RII 


Indicating  concentrations  (at  the  boundary  =?nd  at  the 
interior  respectively) . 

shortest  distance  from  fiber  to  fiber  surface 

diffusion  coefficient 

diffusion  coefficient  of  the  resin 

diffusion  coefficient  of  the  fiber 

average  diffusion  coefficient 

effective  diffusion  coefficient 

j times  the  thickness  of  the  composite  plate 

flux  of  diffusant 

distribution  coefficient 

thermal  conductivity  of  the  matrix 

thermal  conductivity  of  the  fiber 

electrical  conductivity  of  the  matrix 

thickness  of  the  composite  plate 

flux  (general) 

amount  absorbed  at  time  t 

amount  absorbed  after  equilibrium 
has  been  reached 

fiber  radius 

relative  humidity 
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II 
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II 

X 

X 

XX 
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time 

parallel 

perpendicular 

Subscripts 
Symbols  Indicating 

effective 

fiber 

hexagonal 

matrix 

resin 

tetragonal 

parallel 

perpendicular 

first  derivative  of  the  variable  x 
second  derivative  of  the  variable  x 
first  deriva^ve  of  the  variable  y 
second  derivative  of  the  variable  y 
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